arXiv:1504.03690v3 [astro-ph.SR] 15 May 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 18 May 2015 (MN style file v2.2) 

Inefficient star formation through turbulence, magnetic 
fields and feedback 

Christoph Federrath^* 

^Research School of Astronomy and Astrophysics, The Australian National University, Canberra, ACT 2611, Australia 


18 May 2015 


ABSTRACT 

Star formation is inefficient. Only a few percent of the available gas in molecular clouds 
forms stars, leading to the observed low star formation rate (SFR). The same holds 
when averaged over many molecular clouds, such that the SFR of whole galaxies is 
again surprisingly low. Indeed, considering the low temperatures, molecular clouds 
should be highly gravitationally unstable and collapse on their global mean freefall 
time-scale. And yet, they are observed to live about 10-100 times longer, i.e., the SFR 
per freefall time (SFRff) is only a few percent. Thus, other physical mechanisms must 
counteract the quick global collapse. Turbulence, magnetic fields and stellar feedback 
have been proposed as regulating agents, but it is still unclear which of these processes 
is the most important and what their relative contributions are. Here we run high- 
resolution simulations including gravity, turbulence, magnetic fields, and jet/outflow 
feedback. We confirm that clouds collapse on a mean freefall time, if only gravity is 
considered, producing stars at an unrealistic rate. In contrast, if turbulence, magnetic 
helds, and feedback are included step-by-step, the SFR is reduced by a factor of 2-3 
with each additional physical ingredient. When they all act in concert, we find a con¬ 
stant SFRff = 0.04, currently the closest match to observations, but still about a factor 
of 2-4 higher than the average. A detailed comparison with other simulations and with 
observations leads us to conclude that only models with turbulence producing large 
virial parameters, and including magnetic fields and feedback can produce realistic 
SFRs. 

Key words: MHD — ISM: clouds — ISM: kinematics and dynamics — ISM: jets 
and outflows — stars: formation — turbulence 


1 INTRODUCTION 

Why is star formation so inefficient? For example, our en¬ 
tire home galaxy—the Milky Way— only produces about 

2 Mq yr“^ jchoiniuk fc Povichll20li] l. while there is plenty 
of gas available to f orm h undreds of stars per yea r 
lIZudcerman fc Palmer! Il974l : IZuckerman fc Evand Il974ll . 
Based on the range of typical molecular cloud masses 
of ~ 10^-10^ Mq and sizes of ~ 1-100 pc, the average 
cloud mass density p is a few times 10“^^-10~^® gcm~®, 
corresponding to an average particle number density 
of about 10^-10® cm~® for standard molecular compo¬ 
sition. This would render the clouds highly unstable 
and immediately leads to a very short freefall time, 
ts = y^3n/{32Gp) = 0.1-1 Myr, which is at least one order 
of magnitude shorter than the typical lifetime of molecular 
clouds. Thus, other mechanisms than thermal pressure must 
be at work to prevent molecular clouds from collapsing glob- 
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ally and from forming stars at a 10-100 times higher rate 
than observed. 

For a long time, the standard picture of star formation 
was that magnetic fields provide the primary source of sta¬ 
bilizing pressure and tension, which acts against the quick 
global collapse. Only after the neutral species had slowly 
diffused through the charged particles over an ‘ambipolar- 
diffusion’ time-scale of about 10 Myr, would star forma- 
tion proceed in the central regions of magnetized clouds 
llMestel fc SpitzeJll956l : [Mouschovi'^ll976l : IShulll983f . This 
standard theory requires that clouds start their lives dom¬ 
inated by a strong magnetic field, rendering them ini¬ 
tially subcritical to collapse. After about one ambipolar 
diffusion time-scale, magnetic flux was left behind in the 
cloud envelope, while the mass has increased in the cloud 
core such that some stars could form in centre. Thus, 
star formation regulated by ambipolar diffusion predicts 
a higher mass-to-flux ratio in the cores than in the en¬ 
velopes of the clouds, which i s —however—typically not ob- 
served dCrutcher et al.l l2009l : iMouschovias fc Tassisl l2009l : 
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Lunttila et alJl2009l: 

Santos-Lima et al.ll2010l:lLazarian et alJ 

2 OI 2 I: Bertram et al. 

201211 . 


An alternative scenario is that clouds are in fact col¬ 
lapsing globally, thus initially forming stars at a very high 
rate, but that the ionization feedback from massive stars 
eventu ally terminates the global collapse and disperses the 
clouds (IVazauez-SemadeMl2015ll . 


A third alternative is that the observed super- 


sonic 

random motions 

lIZuckerman & Palmer 1974 

Zuckerman & Evand 19741: 

LarsonI Il98ll: ISolomon et al.l 

19871: 

iFalgarone et al.l 

I 992 I: lOssenkoof & Mac Lowl 

2 OO 2 I: 

Hever & BruntI l2004 ISchneider et al. 201ll: 

Roman-Duval et al.l l2011^ 

regulate star formation. In 


this modern picture, turbulence plays a crucial dual role. 
On the one hand, the turbulent kinetic energy stabilizes 
the clouds on large scales and prevents global collapse, on 
the other hand, it induces local com pressions in shocks, 
because the turbulence is supe r sonic dMac Low fc KlessenI 
I 2 OO 4 I : lElmegreen fc Scai3 l2004l : iMcKee fc Ostrikerl l2007ll . 
This generates the initial conditions for star formation, 
because the local compressions typically produce hla- 
ments and dense cor es at the intersections of filaments 
llSchneider et al.ll2012h . 


So which scenario is right or wrong? As it turns out, 
neither of the three is entirely wrong, nor do they explain 
everything by themselves. In the most recent years, we have 
come to a deeper understanding of the relative importance 
of magnetic fields and turbulence through numerical simula¬ 
tions. While the role of turbulence cannot be denied, because 
it naturally explains most of the observed velocity and den¬ 
sity structure of molecular clouds, the role of magnetic fields 
remai ned le ss clear. Simulation work by e .g., IPrice fc Bate! 
ll2007h and iHennebelle fc Tevssieil (l2008l ~) suggested that 
magnetic helds can suppress fragmentation on small scales 
and hence influen ce the shape of the initial mass function 
of stars. Moreover, Wan^_et aljj [201Cll l. |Padoan fc Nordlundl 
ll201lll , iFederrath fc KlessenI 1 2012l ~l and Mvers et al.l ( 2014l l 
showed that magnetic fields reduce the SFR by a factor of 
2-3 compared to the non-magnetized case. Importantly how¬ 
ever, magnetic fields also launch fast, powerful, mass-loaded 
jets and outflows from the protostellar accretion disc. These 
jets and outflows can drive turbulence and alter the cloud 
structure and dynamics so profoundly that the SFR and 
the initial mass function might be even more affected by 
this jet/outflow feedb ack mechanism (|Krumholz et al.ll2014l : 
IFederrath et al.ll2014h . 


The aim of this study is to determine the physical 
processes that make star formation inefhcient. We measure 
which physical mechanisms are relevant for bringing the SFR 
in agreement with observations and we determine their rela¬ 
tive importance, i.e., by what amount they reduce the SFR 
individually and when acting all in concert. 

Section [ 2 ] summarizes our numerical methods and sim¬ 
ulations from which we measure the SFR. Our results are 
presented in Section (3] where we find that purely self- 
gravitating molecular clouds are indeed highly unstable, 
while the step-by-step inclusion of turbulence, magnetic 
fields, and jet/outflow feedback brings the SFR down by one 
to two orders of magnitude when all of these physical mech¬ 
anisms act together. In Section [fl] we review the existing 
literature and compare our simulations and those of other 


groups with observational data for the SFR. We summarize 
our findings and conclusions in Section [S] 


2 NUMERICAL SIMULATION TECHNIQUES 

We use the multi - physi cs, adaptive me sh refinement (AMR 
Berger fc Colell^ 1989l l code FLASH ijFrvxell et al.l l200fll : 
Dubev et al. l20^ i in its latest version (v4), to solve the 


compressible magnetohydrodynamical (MHD) equations on 
three-dimensional (3D) periodic grids of fixed side length L, 
including turbulence, magnetic fields, self-gravity and out- 
flow feedback. The p ositive-definite HLL5R Riemann solver 
ijWaagan et al.ll201ll ') is used to guarantee stability and ac¬ 
curacy of the numerical solution of the MHD equations. 


2.1 Turbulence driving 

Turb u lence is a key for star form a tion (TMac Low fc KlessenI 
2 OO 4 I: Elmegreen fc ScaI3 20oJ: McKee fc Ostrikerl l2007l: 


Krumholz 


Padoan et akll^ldh . so most of our sim¬ 


ulations include a turbulence driving module that produces 
turbulence similar to what is observed in real molecular 
clouds, i.e., drivin g on the largest scales l|Hever et al ] l200fil: 
iBrunt et al.l l2009li and with a power spectrum, S(k) ~ 
co n sisten t with superson i c, compressible turbulence 
jLjmonl Il98ll : iHever fc BruntI l2004l : iRomamDuvaJet^^jJ 
201l|) and confirmed by simulations (IKritsuk et al.l 200?! : 
Federrath et al.l l2010bl : (Pederra^ l2013al 'l. We drive turbu¬ 


lence by applying a sto c hastic Ornstein-Uhle nbeck process 
jEswaran fc Pope 1 19881 : ISchmidt et al.l I 2 OO 6 I I to construct 
an acceleration field Fstir, which serves as a momentum and 
energy source term in the momentum equation of MHD. 
As suggested by observations, Fstir only contains large-scale 
modes, 1 < |k|L/27r < 3, where most of the power is in¬ 
jected at the kinj = 2 mode in Fourier space, i.e., on half of 
the box size. The turbulence on smaller scales is not directly 
affected by the driving and develops self-consistently. The 
turbulence forcing module used here excites a natural mix¬ 
ture of solenoidal and compressible mode s, corresponding 
to a t urbulent driving parameter b = 0.4 l|Federrath et al.l 
l2010bl l. although some cloud-to-cloud variations in this pa¬ 
rameter from b ^ 1/3 (purely solenoidal driving) to 6 ~ 1 
purely compress i ve driving) are expecte d for real clouds 
Price et al.ll201ll : iKainulainen et ^120131 1. 


2.2 Sink particles 

In order to measur e the SFR, we use th e sink particle 
method developed bv iFederrath et al.l l|2010a^ . Sink particles 
form dynamically in our simulations when a local region in 
the simulation domain undergoes gravitational collapse and 
forms stars. This is technically achieved by first flagging each 
computational cell that exceeds the Jeans resolution density, 

(i) 

with the sound speed Cs, the gravitational constant G and 
the local Jeans length Aj. Thus, the sink particle accretion 
radius is given by rgink = Aj/2 and set to 2.5 grid cell lengths 
in order to capture star formation and to avoid artificial 
fragmentation on the highest level of AMR llTruelove et al.l 
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Table 1. Key simulation parameters. 


Model Name 

(1) 

Turbulence 

(2) 

cTt, (km/s) 

(3) 

M 

(4) 

B(pG) 

(5) 

0 

(6) 

Ma 

(7) 

Jets/Out flows 
(8) 

(9) 

SFR(Mo/yr) 

(10) 

SsFR (Mq /pc^ /Myr) 
(11) 

SFRfi 

(12) 

G 

None 

0 

0 

0 

00 

00 

No 

10243 

1.6x10-4 

39 

0.47 

GvsT 

Mix 

1.0 

5.0 

0 

00 

00 

No 

10243 

8.3x10-3 

21 

0.25 

GvsTM 

Mix 

1.0 

5.0 

10 

0.33 

2.0 

No 

10243 

2.8x10-3 

6.9 

0.083 

GvsTMJ 

Mix 

1.0 

5.0 

10 

0.33 

2.0 

Yes 

10243 

1.4x10-3 

3.4 

0.041 

GvsTMJ512 

Mix 

1.0 

5.0 

10 

0.33 

2.0 

Yes 

5123 

1.3x10-3 

3.2 

0.039 

GvsTMJ256 

Mix 

1.0 

5.0 

10 

0.33 

2.0 

Yes 

2563 

8.9x10-3 

2.2 

0.027 


Notes. Column 1: simulation name. Columns 2—4: the type of turbulence driving llFederrath et al1l2010bh . turbulent velocity dispersion, 
and turbulent rms sonic Mach number. Columns 5-7: magnetic field strength, the ratio of thermal to magnetic pressure (plasma /3), and 
the Alfven Mach number. Column 8: whether jet and outflow feedback was included or not. Column 9: maximum effective grid resolution 
(note that refinement is based on the Jeans length with a minimum of 32 cells per Jeans length). Columns 10-12: the absolute SFR, 
the SFR column density, and the SFR per mean global freefall time. The standard four simulation models are in the first four rows. 
Two additional simulations with lower resolution, but otherwise identical to GvsTMJ, are listed in the last two rows, in order to check 
convergence of our results for the SFR. 


I 1997 I ). If the gas density in a cell exceeds paink, a spherical 
control volume with radius rgink is constructed around that 
cell and it is checked that all the gas within the control vol¬ 
ume is Jeans-unstable, gravitationally bound and collapsing 
towards the central cell. A sink particle is only formed in the 
central cell of the control volume, if all of these checks are 
passed. This avoids spurious formation of sink particles and 
guarantees that only bo und and collapsing gas forms stars 
jFederrath et al.ll2010ah . which is important for accurately 
measuring the SFR. 

On all the lower levels of AMR (except the highest level, 
where sink particles form), we use an adaptive grid refine¬ 
ment criterion based on the local Jeans length, such that 
Aj is always resolved with at least 32 grid cell lengths in 
each of the three spatial directions of our 3D domain. This 
resolution criterion is very conservative and computation¬ 
ally costly, but guarantees that we resolve turbulence on the 
Jeans scale (iFederrath et al.ll2011^ . potential dyna mo ampli¬ 
ficati on of the magnetic field in the dense cores (ISur et al.l 
I2OIOI) . and capture the basic structure of accretion discs 
forming on the smallest scales (IFederrath et al.ll2014l) . If a 
cell within the accretion radius of an existing sink particle 
exceeds Psink during the further evolution, is bound to the 
sink particle and is moving towards it, then we accrete the 
excess mass above Psink on to the sink particle, conserving 
mass, momentum and angular momentum. We compute all 
contributions to the gravitational interactions betw een the 
gas o n the grid (with the iterative multigrid solver bv I Ricked 
I2OO8I) and the sink particles (by direct summation over all 
sink particles and grid cells). A second-order leapfrog inte¬ 
grator is used to advance the sink particles on a timestep 
that allows us to resolve c lose and highly eccent ric orbits 
(for details, see the tests in iFederrath et aDl2010ah . 

2.3 Outflow/Jet feedback 

Powerful jets and outflows are launched from the protostellar 
accretion discs around newborn stars. These outflows carry 
enough mass, linear and angular momentum to transform 
the structure of their parent molecular cloud and to po¬ 
tentially control star formation itself. In o rder to take this 
most important mechanical feedback effect (iKrumholz et al.l 


I 2 OI 4 I) into account, we recently extended the sink particle 
approach such that sink particles can launch fast collimated 
jets together with a wide-angle, lower-speed outflow compo¬ 
nent, to reproduce the global features of observed jets and 
outflows, as well as to be consistent with high-resolution 
simulations of t he jet launching proc ess and with theoreti¬ 
cal predictions (IFederrath et al.ll20l4l . The most important 
feature of our jet/outflow feedback model is that it converges 
and produces the large-scale effects of jets and outflows al¬ 
ready with relatively low resolution, such as with sink par¬ 
ticle radii rsink ~ 1000 AU used in our star cluster simula¬ 
tions here. Our module has been carefully tested and com¬ 
pared to previous implementations of jet/outflow feed back 
such as the models imple mented in IWang et al.l (I2OIOI) and 
ICunningham et al.l (1201 j) . The most important difference to 
any previous implementation is that our feedback model in¬ 
cludes angular momentum transfer, reproduces the fast col¬ 
limated jet component and dem onstrated convergence (for 
details, see IFederrath et ^I20l4l . 

2.4 Simulation parameters 

All our simulations share the same global properties: a 
cloud size L — 2 pc, a total cloud mass M — 388 Mq and 
a mean density po = 3.28 x 10“^'^gcm“®, resulting in a 
global mean freefall time tfr = 1.16 Myr. Models including 
turbulence have a velocity dispersion (j„ = lkms“^ and 
an rms Mach number of AI = 5, given the sound speed 
Ca = 0.2kms“^, appropriate for molecular gas with temper¬ 
ature T — lOK over the wide range of dens ities that lead 
to de nse core and eventually star formation (lOmukai et al.l 
I2OO5I) . Finally, models including a magnetic field start with 
a uniform initial field of B = 10 pG, which is subse¬ 
quently compressed, tangled and twisted by the turbu¬ 
lence, similar to how it would be structured in real molec¬ 
ular clouds (IPlanck Collaboration et al.l I2OI5I ). The mag¬ 
netic field strength, the turbulent velocity dispersion and 
the mean density all follow typical values derived from 
observations of clouds with the given phy sical properties 
(iFalgarone et al.l 1 19921 : ICrutcher et al.l 120101 ) . This leads to 
the dimensionless virial ratio Ovir = 1.0 (also typical for 
molecular clouds in the Milky Way; see iKauffmann et aP 
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Figure 1. Column density projections at the end of each simulation model: Gravity only (top left), Gravity vs. Turbulence (top right), 
Gravity vs. Turbulence + Magnetic Fields (bottom left), and Gravity vs. Turbulence + Magnetic Fields + Jet/Outflow Feedback (bottom 
right). The time to reach a realistic SFE = 20% with small star clusters having formed, is shown in the top right corner of each panel. 
The time increases significantly for each model that adds physical processes opposing gravitational collapse. This demonstrates that only 
the combination of turbulence, magnetic fields and feedback can produces realistic SFRs, which is quantified below and summarized in 
Table [T] An animation of these still frames is available in the online version. 


l2013h and to a plasma beta parameter (ratio of thermal to 
magnetic p ressure) 0 = 0.33 or a n Alfven Mach nnmber 
Ma = 2.0. iFalgarone et al.l (l2008l ) find an average Alfven 
Mach number of about Ma = 1-5 in 14 different star¬ 
forming regions in the Milky Way. Thus, the assumed mag¬ 
netic field in our simulation models is very close to the values 
typically observed in molecular clouds and cloud cores. 

We run four basic models, which—step by step—include 
more physics. In the first simulation we only include self¬ 
gravity with a given initial density distribution resembling 
molecular cloud structure, but we do not include any tur¬ 
bulent velocities or magnetic field. In the second model, we 
include a typical level and mixture of molecular cloud tur¬ 
bulence (see mi- The third model is identical to the sec¬ 


ond model, but adds a standard magnetic field for the given 
cloud size and mass. Finally, the fourth model is identical 
to the third model, but additionally includes jet and out¬ 
flow feedback (see 1 12.311 . These four basic models were all 
run with a maximum effective grid resolution of 1024® cells. 
Their key parameters are listed in Table [T] We also run two 
additional models, which are identical to the fourth model 
(with jet/outfiow feedback), but have a lower maximum ef¬ 
fective resolution of 512® and 256® cells, respectively, in or¬ 
der to check numerical convergence of our results for the 
SFR. Those models are listed in the bottom two rows of 
Table [T] 
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3 RESULTS 


3.1 Cloud structure and stellar distribution 

Figure [T] shows column density projections of our four basic 
models from Table [T] The simulation that only includes self¬ 
gravity (no turbulence, no magnetic fields and no feedback) 
is shown in the top left-hand panel. The gas structures we 
see in the figure resemble the typical distribution found in 
gravity-only simulations such as simulations of pure dark 
matter in the cosmological case. Gas from the voids con¬ 
tinuously falls on to dense filaments. Stars form inside the 
filaments and in particular where filaments intersect. This 
seems to be true also fo r real molecular clouds to some degree 
l|Schneider et al.ll2012h . but the inter-filament gas is much 
more disturbed and turbulent in real molecular clouds com¬ 
pared to the gravity-only case. We stop the simulation when 
20 % of the gas is wound up in stars, i.e., the star formation 
efficiency SFE = Mstars/(Mstars -b Mgas) = 20%, which is 
already reached within a fraction of a global mean freefall 
time, t = 0.61 tff, for the gravity-only simulation. Thus, the 
whole cloud forms stars in roughly a global freefall time, 
which is the expected—and unrealistic—outcome if no other 
physics but gravity is considered. 

When we add a realistic level of turbulence to the cloud, 
shown in the top right-hand panel of Figure [U we find that 
the time to reach SFE = 20% increases to t — 1.31 fg or 
1.5 Myr. The overall structure is now much closer to what 
we see in real molecular clouds, especially the turbulent ve¬ 
locity dispersion measured in observations is approximately 
reproduced in this model. However, the SFR is still an or¬ 
der of magnitude too large compared to typical observations 
and there is no magnetic field present, contrary to what is 
observed JCrutcher et akllioiol i. 

Thus, we add magnetic fields in the simulation shown 
in the bottom left-hand panel of Figure [T] and see that 
the end time is shifted by slightly more than a fac¬ 
tor of 2, f = 2.91 fff = 3.4 Myr, compared to non- 
MHD case. Magnetic fields thus reduce the SFR by 
about a factor of 2-3, as found in previous simulations 
jlWaneet_^ 20101: ^doan fc Nordlundll201ll : |Padoan et al.l 


l2012l : lFederrath fc KlesserJl2012l i. Although we added a con¬ 

siderable set of physics (gravity, turbulence and magnetic 
fields) in this simulation model, the SFR is still much 
higher than inferred from observations, but many spatial 
and morphological features seen in real star-forming molec¬ 
ular clouds are well reproduced with this set of physics. 

Finally, we add jet and outflow feedback in the simu¬ 
lation shown in the bottom right-hand panel of Figure [T] 
The overall cloud structure is similar to the turbulent MHD 
case without feedback, but the SFR is further reduced by 
about another factor of 2. Now it takes 6.5 fg = 5.6 Myr to 
reach a reasonable SFE = 20% in this small cluster-forming 
region of a molecular cloud. The most interesting morpho¬ 
logical difference compared to all the no-feedback cases is 
that the oldest stars have wandered away from their for¬ 
mation sites and are now distributed over a much wider 
volume. Newborn stars are still always located on or very 
close to high-density peaks. A total of 37 stars have formed 
in the whole field, while some of those stars might actu¬ 
ally represent close binaries or small multiple systems given 
the limited resolution available in these calculations. While 
the absolute fragmentation is not fully converged in these 
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Time per freefall fime (f/f,,) 


Figure 2. Time evolution of the SFE (top panel) and the SFR per 
freefall time, SFRg = d(SFE)/ci(t/tg) (bottom panel) for each of 
our four standard models from Table [T] We immediately see the 
drastic reduction in the SFR when turbulence, magnetic fields 
and jet/outfiow feedback together counteract gravity, in which 
case the SFR is reduced to values consistent with the observed 
range of SFRs, shown as the hatched region. 


simulations, we emphasize that the SFR and the amount 
of mechanical feedback is well captured and is converged 
with our numerical resolution. This is shown in Table [T] 
where in the last three columns, we compare this feed¬ 
back model (GvsTMJ) with two lower-resolution equivalents 
(GvsTMJ512 and GvsTMJ256), demonstrating convergence 
of the SFR to within 5%. 


3.2 Star formation rate 

Now we concentrate on the quantitative determination of 
the SFR. Figure [2] shows the time evolution of the SFE 
(top panel) and the SFR (bottom panel) in our four 
fiducial simulation models. While the gravity-only model 
(dot-dashed line) turns all the gas into stars in about a 
global mean freefall time, the S FR per freefall time defined 
as SFRg = d(SFE)/d(f / fg) dKrumholz fc McKed l2005l: 
iFederrath fc KlessenI [2012 ) is reduced by about a factor of 
2 when turbulence is included (dashed line). The SFRg is 
reduced by at least another factor of 2 when turbulence and 
magnetic fields are present (dotted line), and by another 


© 0000 RAS, MNRAS 000, 000-000 




































6 Federrath 


factor of 2 when jet and outflow feedback is included (solid 
line). 

The hatched region in the bottom panel of Figure [2] 
shows the observed ra nge of SFRg based on t he observa¬ 
tional data compiled in iKrumholz fc TanI ll2007ll . Our feed¬ 
back model eventually settles into that observed range. 

A very interesting and noteworthy feature of our feed¬ 
back model is the evolution of the SFR with time, in par¬ 
ticular the fact that the SFR slowly but steadily increases 
until about t ~ 4.0 Myr and then drops significantly by al¬ 
most an order of magnitude until t ~ 5.3 Myr, followed by 
a period of nearly constant SFR. This is a manifestation 
of self-regulation by feedback. As more material is accreted 
and the SFR increases, the amount of gas being re-injected 
into the interstellar medium in the form of fast mass-loaded 
jets and outflows increases proportionally. Eventually, the 
feedback enhances the turbulence and the amount of kinetic 
energy in the system, thereby increasing the virial parame¬ 
ter, CTvir = 2i?kin/Egrav, such that Star formation is rapidly 
quenched. This quenching of the SFR in turn reduces the 
amount of feedback such that the system eventually set¬ 
tles into a self-regulated state of star formation in which 
any intermittent increase in accretion triggers feedback that 
regulates the SFR down to a nearly constant level. 

The last three columns in Table[T]list the time-averaged 
SFRs in each simulation model. For the gravity-only model, 
we find unrealistically high SFRff = 0.47, while the feed¬ 
back model including turbulence and magnetic fields has a 
time-averaged SFRfj = 0.041. Thus, the combination of tur¬ 
bulence, magnetic fields, and jet/outfiow feedback reduces 
the SFR by more than an order of magnitude and brings 
the SFR into the observed range. 


4 COMPARISON WITH OBSERVATIONS 


In Figure [3] we show how our simulations and those by 
other groups compare to observations of the SFR. The 
observational data are a collection of star-forming Milky 
Way clouds, molecular clumps and y o ung stellar o bjects 
(YS Os) from Heiderman et al.l (l201oh . iLada et ^ (l2010l l 
and iGutermuth et al.l (l2 01lll. as well as the Cen tral Molecu¬ 
lar Zone (CMZ) from lYusef-Zadeh et al.l (l2009tl . The extra- 
galactic data of disc (D) and starb urst (SB) galaxie s at low 
and h igh redshift were compiled in iKrumholz et ^ (l2012al . 
l2013h . and summariz ed and extended b y the Small Magel¬ 
lanic Cloud (SMC) in iFederrathI ll2013bf l. 

The top panel of Figure [5] shows the classical Kennicutt- 
Schmidt diagram, i.e., a plot of the SFR surface density 
(Esfr) versus the gas surface density (Egas). The bottom 
panel shows a more r ecent representa t ion of a star forma¬ 
tion law introduced in iKrumholz et ^ (l2012ah . where Esfr 
is plotted as a function of Egas/tfr, i.e., Egas divided by the 
global mean freefall time tfr, which exhibits a better cor¬ 
relation than the classical Kennicutt-Schmidt relation. The 
remaining scatter in this ne w relation is discussed and ex¬ 
plained in IFederrathI ll2013bf ). by variations in the turbulent 
properties—primarily in the turbulent Mach number—of the 
molecular clouds forming stars in the observations. 

In both panels of Figure [3] we superimpose simulation 
data on to the o bservational dat a points . The crosses are 
simulations from iKrumholz et al.l (l2012bl l with the highest 


to the lowest Esfr produced in their simulations without 
turbulence and isolated gravity boundaries (top cross), in¬ 
cluding turbulence and periodic gravity boundaries (middle 
cross), and finally adding outflow feedback (bo t tom c ross). 
The circles show simulations from iMvers et aO ll2014h . The 
top circle is a purely hydrodynamical simulation and the 
bottom circle included a magnetic field, yielding a re¬ 


tent with IWang et al. 

(l2010fl. iPadoan & Nordlundl 1 I 2 OIIII. 

IPadoan et al.l (20121. 

Federrath & KlessenI (20121 and this 

work. Both the Krumholz et al.l (2012bll and Mvers et al.l 


density and extre mely high Egas/ I ff. Th e squares in Figure [3] 
are simulations bv IPrice fc Batel (l2009l l with a strong mag¬ 
netic field. The square with slightly lower Esfr is from their 
simulation model that included radiative heating, while the 
other one did not. From this, we can conclude that radia¬ 
tive heating—at least from low mass stars—only reduces the 
SFR by < 20%. 


The star symbol s in F igure B] show a series of simula¬ 
tions by IWang et al.l 1 I 2 OIOI I. similar to ours, but at higher 
mean column density. The stars from high to low Esfr 
correspond to gravity-only, plus turbulence, plus magnetic 
fields and plus outflow feedback, respectively. Our set of ba¬ 
sic models is shown as diamonds with error bars based on 
the time-vari ations of the SFR s hown in Figured] When we 
compare the IWang et al.l 1 I 2 OIOI I series of simulation mod¬ 
els with the ones in this work, we find the same basic 
trend. Adding turbulence, magnetic fields and outflow feed¬ 
back reduce the S FR step by step. Ho wever, the minimum 
SFRff obtained in I Wang et al.l (I 2 OI 0 II was SFRff = 0.10, 
which is significantly larger than our best model, which has 
SFRff = 0.0411 q q 2 i. Given the higher virial parameter of 
Ovir = 2£'kin/Egrav = 1.5 Compared to our value (cTvir = 1) 
and their very strong magnetic field (lOO/rC compared to 
our 10 /iC), we would have expected that the IWang et al.l 
1 I 2 OIOII simulation would form stars at a lower rate, because 
of the stronger support from the turbulence and the stronger 
support from the magnetic field than in our case. The most 
likely reason for this discrepancy is that IWang et al.l (I 2 OI 0 II 
use a centrally concentrated density profile with p ~ r~^ 
as an initial condition such that star formation in the cen- 
tre of the cloud is mu ch more violent than in the outskirts 
llCirichidis et al.ll201ll '). 

A colliding flow simulation by iBaneriee et al.l (l2009l j 
with magnetic field is also shown in Figure [3] as an arrow. 
In order to compute Esfr, Egas and ts for their simula¬ 
tion mode l, we used the a v erage SFR = 5 x 10~^ Mq yr“^ 
quoted in IBaneriee et al.l (l2009l l. the average cloud mass 
of 1.5 X 10“* Mq, the average cloud area of 7r(45pc)^ and 
the average cloud thickness of 2 pc. This leads to Esfr = 
0.079 Mq pc“^ Myr“^, Egas = 2.4 Mqpc“^, tg = 7.5 Myr, 
and thus SFRff = 0.25. 

Finally, the triang le in Figure [3] shows the simulation 
from IClark et al.l (I2OO5I I. which reaches low SFRs despite the 
fact that neither magnetic fields nor feedback were included, 
but just turbulence. The explanation f or this beha v iour i s 
that the virial parameter is Qvir = 4 in IClark et al.l 1I2OO5II . 
Thus, their Qvir is very high, such that large fractions of the 
cloud are essentially unbound and cannot form stars. This is 
in line with the theoretical predic tion that the SFR drops ex- 
ponentially with increasing llKrumholz fc McKedlioOsI : 
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Figure 3. SFR column density (Ssfr) ^ function of the gas surface density, Egas (classica l Kennicutt-Schmidt diagram, top panel) 
and SsFR as a function of the gas surface density divided by the mean freefall time, Egas/^ff llKrumholz et al]|2012al : lFederrathll2Q13bl . 
bottom panel). The gre y data points are fr om observations in the Milky Way, the SMC, and in disc and starburst galaxies at low and high 
redshift summarized i n l Federrathl 1201^). The coloured and labe lled symbols are fr om nu merical simulations : iKrumho lz et al!| l|2012bl. 
crosse s), I Mvers et al.l i|2Q14l . circles) JPrice &: Batd l|2009l . squares). IWang et al.l l|2010l . stars). Iciark et"^ l|20Q5l . triangle). iBaneriee et aP 
<l2009l . arrow), and from this work (diamonds). The simulations use a wide range of different setups, initial conditions, physics and 
numerical methods. The general trend is that most simulations produce SFRs about one to two orders of magnitude too high compared 
to the typical values found in the observations. Only simulations that include turbulence, magnetic fields and jet/outflow feedback or 
have very strong turbulence (producing very large virial ratios Ovir = 2i?kin/®grav -J 4) are consistent with the observations. 


IPadoan et al.l LoiSl : iFederrath fc KlesserJ I2OI2II . However, 
only few clouds have such high virial p arameters and the 
SFR per freefall time is still quite high in IClark et ahl (120051 1 
with SFRff ~ 0.10, despite the large avu- Thus, magnetic 
fields and feedback seem to be essential physical ingredients 
to reach realistic SFRs. 


5 SUMMARY AND CONCLUSIONS 

We presented high-resolution hydrodynamical simulations of 
star cluster formation including turbulence, magnetic fields 
and jet/outflow feedback. Although our simulation models 
with turbulence, magnetic fields and jet/outflow feedback 
produce SFRff O.Od/lQ Qj, which is currently the clos¬ 
est available match to observations, they still overproduce 
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stars by a factor of 2-4 compared to the observed average 
SFRfj 0.01. This discrepancy might be resolved by con¬ 
sidering other types of feedback in addition to jet/outflow 
feedback, such as radiation pressure, which seems to be ca- 
pable of reducing the S FR further by another factor of 2 
llMacLachlan et aDl2015l . note however that radiative heat¬ 
ing from low- and intermediate-mass stars does not affect 
the SFR significantly; cf. Section 0. Our results also indi¬ 
cate that the typical level of turbulence and magnetic fields 
in molecular clouds may be higher than previously thought. 
We conclude that only the combination of strong turbulence 
with virial parameters above unity, strong magnetic fields 
and mechanical plus radiative feedback can produce realis¬ 
tic SFRs consistent with observations. 
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